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Abstract 

We study spinodal phase separation in unstable thin liquid films on chemically disordered 
substrates via simulations of the thin-fihn equation. The disorder is characterized by immobile 
patches of varying size and Hamaker constant. The effect of disorder is pronounced in the 
early stages (amplification of fluctuations), remains during the intermediate stages and van- 
ishes in the late stages (domain growth). These findings are in contrast to the well-known 
effects of quenched disorder in usual phase-separation processes, viz., the early stages remain 
undisturbed and domain growth is slowed down in the asymptotic regime. We also address the 
inverse problem of estimating disorder by thin-film experiments. 
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Introduction 



Consider a system which is rendered thermodynamically unstable by a sudden change of param- 
eters, e.g., temperature, pressure. The subsequent evolution of the system is characterized by the 
emergence and growth of domains enriched in the preferred phases.- This domain growth process 
is of great importance in science and technology. In this context, an important paradigm is the 
phase-separation kinetics of an initially homogeneous binary (AB) mixture. If the initial concen- 
tration fluctuations grow spontaneously, the evolution is referred to as spinodal decomposition. 
The segregating mixture evolves into coexisting domains of A-rich and B-rich phases. In pure 
(disorder- free) systems, these domains are characterized by an increasing length-scale, L{t), which 
grows as a power-law, L{t) ~ t'^ in time t. The exponent depends on the transport mechanism. 
For diffusive transport, L{t) ~ f^/-^, which is known as the Lifshitz-Slyozov (LS) growth law and 
is based on the evaporation-condensation mechanism.- 

In recent work,- we have emphasized the analogies and differences between usual phase- 
separation kinetics and the spinodal phase separation (SPS) dynamics of an unstable thin liquid 
film (< 100 nm) on a pure substrate. Typically, random fluctuations in the free surface of initially 
flat films grow and evolve into distinct morphological phases, viz., a thinner low-curvature flat film 
phase, and a thicker high-curvature droplet phase. Spinodal growth occurs when the excess inter- 
molecular free energy AG(/j) shows a minimum and the spinodal parameter, (9^AG/ dh^\h=hQ < 0. 
Here, h is the film thickness, and ho is the average thickness. In |[figure][l][]l[ we show a 
typical form of AG vs. h and the corresponding spinodal parameter, d^AG/ dh^ vs h. The double- 



tangent construction for AG in [figure][l][] 1 shows that the film segregates into phases with h = hm 
and h = oo. (This should be compared with phase- separation problems, which are described by a 
double- well potential, and there are two possible values for the equilibrium composition. ^•^) In an 
alternative scenario, commonly known as true dewetting, the thin film breaks up into dewetted dry 
spots (with h = 0) surrounded by undulating liquid ridges and there is no SPS.^ Our recent work- 
on SPS on homogeneous substrates has shown that the number density of local maxima in the 
film's surface (defects), A^^, is a suitable marker to identify the early, late and intermediate stages 
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Figure 1: Variation of the free energy per unit area (AG, solid line) and spinodal parameter 
(d^Mjjdh^, dashed line) with film thickness. The parameter values are R = —0.1 and 5 = 10 
nm. Spinodal phase separation takes place for h > 12.5 nm, where d^AG/ dh^ < 0. 

of SPS in thin films. We have also established that the kinetics of the early stages can be described 
by a universal evolution of this marker with a suitably rescaled time, and the late stages can be 
described by the LS growth law. 

The above discussion has focused on SPS in pure systems with homogeneous substrates. How- 
ever, real experimental systems never satisfy this condition. The substrates contain physical and 
chemical impurities (disorder), which can be immobile (quenched) or mobile (annealed). We have 
some understanding of the alternative scenario of true dewetting on disordered substrates,- but the 
problem of SPS with disorder is yet to be addressed. This letter investigates SPS in unstable thin 
liquid films on substrates with quenched chemical disorder, and highlights novel features which 
arise due to the presence of disorder. The inverse problem of estimating the disorder through 
thin-film experiments is also addressed, as this would be very useful for experimentalists. 
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Model and Simulation 



The starting point of our study is the thin-film equation which models the evolution of the film's 
surface h{x, t) in supported thin liquid films. This is derived by simplifying the equations of motion 
via the lubrication approximation—. The resulting equation is analogous to the Cahn-Hilliard (CH) 
equation- of phase separation with an /z-dependent mobility M{h) = h? (corresponding 
to Stokes flow with no slip). The total free energy is F^/i] = /[AG(/z) + 7(V/z)2/2] Jjc = Fg + Fi, 
where F^ is the overall excess free energy and Fi is the interfacial free energy. The quantities 7 
and }i refer to surface tension and viscosity of the liquid film, respectively. The corresponding CH 
equation is 



-h{x,t)=V- 



MV 



8Fs 



3ix \ dh ' 



(1) 



where all gradients are taken in the plane of the substrate. The potential AG(/7) usually combines a 
long-range attraction and a short-range repulsion. The results presented here correspond to a long- 
range van der Waals attraction due to the substrate, and a comparatively short-range van der Waals 
repulsion provided by a nano-coating on the substrate: AG = —Ac/llnh^ — As / \2Tt{h -\- 5)^. 



This is the potential shown in |[figure][l][]l[ Here, As and Ac (= RAs) are the effective Hamaker 
constants for the system, which consists of the fluid above the film, the liquid film, and a solid 
substrate (s) or coating material (c). The thickness of the nano-coating is 5. 



We can reformulate [equation] [1][]1 in a dimensionless form as follows: 



—HiXJ)=V. 
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In |[equation] [2] [] 2] H = h/ho, where Hq is the mean film thickness; X =x/^, where ^ = i^Tiyj \As | ) ^'^^ /iq 
is the characteristic scale for the van der Waals case; and T = t/x, where T = {llTf-iiyhyA^^ . The 
excess free-energy term now has the form 



Inhl (9 AG 
I A J dH 
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(3) 
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where D = S/hQis the nondimensional coating thickness. The linear stability analysis of [[equation] [2] []2 
for fluctuations about H = I predicts a dominant spinodal wave of wave- vector kM and time-scale 
Tm with 

kM 





/ nhl d^AG 




-=\\ 
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H=l 
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(4) 



For a homogeneous substrate, the parameters in [[equation] [2] [] 2|[equation] [3] [] 3| are spatially 
uniform. We model the quenched disorder via fixed patches of varying size P with relative 
Hamaker constant R. The mean values of P and R are Pm and Rm- These parameters are uni- 
formly distributed in the intervals [P,,, - bP.Pm -f SP] and [i?,„ - bR.Rm -V 5R], with 5P/P,n = Pd 



and 5R/R,n = Rd- In [[figure] [2] []2[ we plot the spinodal parameter {2nhQ/ \ As \)d AG/dH vs. 
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Figure 2: Variation of dimensionless spinodal parameter with nondimensional film thickness for 
three values of = -0.15,-0.1 and -0.05. 



H for three representative values of R. These correspond to the mean value R,„ = —0.1, and the 
extreme deviations for i?^ = 0.5. Notice that a small change in R results in a major change of kM 
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(or Lm) and Tm in [[equation] [4] []4 

We numerically solve [equation] [2] [] 2 in J = 2,3 starting with a small- amplitude (~ 0.01) 
random perturbation about the mean film thickness H = I. The system size in J = 2 is hLm, 
where Lm is the dominant wavelength for R = {n ranges from 16 to several thousands). The 
system size in J = 3 is (16Lm) ■ We apply periodic boundary conditions. A 512-point grid per 
Lm was found to be sufficient when central-differencing in space (with half-node interpolation) 
was combined with Gear's algorithm for time-marching, which is convenient for stiff equations. 
The parameters D, R,n and R^ were chosen so that the film is spinodally unstable at // = 1 (i.e., 
a > 0) for all values of R. A typical spatial variation of R with R^ = —0. 1 and R^ = 0.5 is shown 



in 



[figure][3][]3 P,n was taken as /Lm, where / varies from 1/16 to 16. Thus, the size of the 



patches with fixed R varied from being much smaller than the spinodal length-scale to being much 
bigger than it. 
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Figure 3: Schematic of a chemically disordered substrate. The disorder is quenched and is charac- 
terized by R,„ = —0. 1 ( horizontal dashed line), Rd = 0.5, Pm = Lm /16 and Pj = 0.5. 
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Results and Discussion 



All the three markers of the process, viz., number density of defects (A'^;), morphology and energy 
are affected by the strength (Rj) as well as the length-scale (Pm) of the chemical disorder. Notice 
that Nj^ characterizes the typical domain size. In [figure] [4] []4 we plot Nj vs. T on a log-log 
scale. The effect of disorder is pronounced in the early stages, persists during the intermediate 
stages, and vanishes in the late stages. The evolution of A^^ briefly follows the results of the pure 
case (R^ = 0), but is seen to diverge while still in the early stages. However, surprisingly, the plot 
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Figure 4: Variation of the number density of defects {N^) with nondimensional time (T) for a 
d = 2 system of size 2048 Lm and D = 0.5. Curves 1 to 6 present results for R^ = —0.1 and 
Rd = 0.50, 0.25, 0.10, 0.05, 0.03 and 0.01, respectively. The LS -1/3 slope for late stages is 
shown by the dashed line. The dotted curve corresponds to the case with zero disorder (base 
case). The inset shows the magnified intermediate stage for Rj = 0.50, 0.05 and 0.01. The other 
parameter values are P,n = Lm /16, = 0.5. 



of N^i ws. T reverts to the pure case in the late stages. This is the central result of this letter, and 
will be discussed shortly. The deviation in the intermediate stages is seen in the magnified view 
in 



[figure] [4] []4 A decrease in the disorder strength delays the divergence, but even a weak disor- 



der amplitude (Rj = 0.01) results in a sizable split (see [figure][4][]4 1. The split occurs because the 
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growth of initial fluctuations is drastically amplified by the presence of disorder ( [[figure] [5] []5 1. 



Consider the growth of initial fluctuations, H = I + 6 {X ,T) , in [equation] [2] [] 2 On Fourier- 
transforming 6{X,T) and integrating the resultant 0{k,T) over a uniform disorder distribution on 
[^m(l ~Rd),Rm{l +Rd)\, we obtain 

e{kj)= exp [k^ { a {R,n) -k^]T] ''''If^^^^^^^^ il 0) , (5) 

where b = 1 — (1 + The factor which represents the effect of the disorder, sinh(cr) / (cT), 

— 7- 1 as — 7- 0, but amplifies the growth of initial fluctuations for nonzero Rd- In turn, this increases 



the sub-structure in the growing profile (see frames at T = 0.052, 10 in [figure] [5] []5 1 - resulting 
in a slower decrease of vs. T. The growing fluctuations are saturated in the late stages when 
domains are formed at the smaller value of the equilibrium thickness < 1 . Of course, the height 
field continues to grow for H > I. These domains coarsen with time and obey the LS growth law. - 
The presence of disorder causes local fluctuations in the domain structure - see frame at T = 10000 
in 



[figure][5][]5 



In the late stages, the growth law is universal in [[figure] [4] []4| (A^^ ~ T ^/^), regardless of the 



presence of disorder. This should be contrasted with the usual phase- separation problem, where 
coexisting domains are trapped by disordered sites.—"— The subsequent coarsening proceeds by 
thermally-activated barrier hopping, leading to an asymptotically logarithmic growth law. In the 
SPS problem discussed here, the defects consist of "interfaces" between flat domains at H = Hi„ 
and growing hills with // — )• oo. These defects do not become trapped as the disorder scale becomes 
irrelevant in comparison to the diverging height of the "interfaces". This late-stage universality has 
important implications for experimentalists. As mentioned earlier, it is not possible to eliminate 
disorder from experimental systems. However, our results show that disorder is irrelevant in the 
late stages, and all experiments will finally recover results for the pure case. The differences 
between the pure and disordered systems only show up in the early stages, as quantified above. 
The evolution of the total free energy is also consistent with the evolution of the defect density 
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Figure 5: Morphological evolution of the film at different nondimensional times, as specified. We 
show results for R = —0.1 and R^i = 0.5,0.05 and the pure case. The other parameter values are 



the same as in [figure] [4] []4 



( [[figure] [6] []6 ). As expected, the overall free energy diminishes with time. The decay in the late 
stages is a universal power-law, which is independent of the disorder amplitude. The difference 
between the pure and disordered cases is seen in the early stages, with faster decay for larger 
disorder values, as the growth of initial fluctuations is speeded up by disorder. 

Next we focus on the effect of varying the patch size (Pm^Pd) on SPS. The early-stage kinetics 
exhibits three qualitatively distinct behaviors depending on the ratio / {Pm = /Lm) ■ Notice that the 
limit Pm — > system size corresponds to the pure case. For small values of / (/ < 1/2), the disorder 



leads to splitting as described earlier [Type A curves in [[figure] [7][]7| ^a)]. For intermediate values 
of/(l/2</<2), the disorder leads to a distinctive staircase behavior for A'^ vs. T [Type B 
curves in 



[figure] [7] []7 a)]. In this case, the patch size is comparable to the spinodal wavelength. 
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Figure 6: Variation of the total free energy with nondimensional time for a system of size 2048 
Lm and D = 0.5. The curves present results for = —0.1 and = 0.50, 0.25, 0.1. The other 
parameter values are = Lyi^/16, = 0.5. The dotted curve corresponds to the case with zero 
disorder. The inset shows the magnified view for energy variation in the early stages. 

resulting in a stick-slip dynamics of the height field. For larger values of / (/ > 2), the disorder has 
little effect and the early-stage kinetics resembles that for homogeneous substrates [Type C curves 
[figure] [7] []7 "a)]. As expected, the late-stage dynamics is again universal with A^^ 



in 



do not show this regime in [[figure] [7][]7| ^a). A similar behavior is also seen in J = 3 simulations 
- in 



[figure] [7] []7 ^b), we plot vs. T for this case. The simulation details are provided in the 
figure caption. 

As the spinodal length-scale can be controlled by changing the mean film thickness Qm °^ h^, 
where Im is the dimensional length-scale), one can use the present findings to address the inverse 
problem of assessing the disorder by thin-film experiments. In thick films, we see Type A behavior 
because / ^ 1 . A reduction in the film thickness such that 1/2 < / < 2 results in Type B behavior. 
A further reduction in thickness (/ > 2) results in Type C behavior. We can use Type C curves to 



first estimate Lm, and then R^ from [equation] [4] []4 Thus, we can estimate the disorder amplitude 



to an order of magnitude. An average value found by monitoring the spinodal length-scale when 
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Figure 7: Variation of the number density of defects with nondimensional time in (a) d = 2 and 
(h) d = 3. The parameters are Rm = —0. 1, = 0.5. The different curves are obtained by varying 
the mean patch size Pm = /Lm- We classify the curves as follows: Type A (/ < 1/2), Type B 



Type B behavior starts and when it ends can fine-tune this estimate. The strength of the disorder 
can be found by any Type C result. 

Conclusions 

In summary, we conclude that quenched chemical disorder has a pronounced effect on the early 
and intermediate stages of morphological phase separation in thin liquid films. However, the late- 
stage kinetics is universal and follows the Lifshitz-Slyozov (LS) growth law. These findings are in 
sharp contrast to the effects of quenched disorder seen in usual phase-separation processes, viz., 
the early stages remain undisturbed and domain growth is slowed down in the asymptotic regime. 
The early-stage kinetics shows different qualitative behavior, depending on the relative sizes of 
disorder patches and the spinodal length-scale. Thus, the inverse problem of estimating disorder 
by thin-film experiments can also be addressed. 



(l/2</<2),andTypeC(/>2). 
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